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Abstract 

In the design and analysis of robust control systems for uncertain plants, the technique of formulating what 
is termed an "M-A model" has become widely accepted and applied in the robust control literature. The ”M" represents 
the transfer function matrix M(s) of the nominal system, and "A" represents an uncertainty matrix acting on M(s). The 
uncertainty can arise from various sources, such as structured uncertainty from parameter variations or multiple 
unstructured uncertainties from unmodeled dynamics and other neglected phenomena. In general, A is a block diagonal 
matrix, and for real parameter variations the diagonal elements are real. As stated in the literature, this structure can 
always be formed for any linear interconnection of inputs, outputs, transfer functions, parameter variations, and 
I^rtuibations. However, very little of the literature addresses methods for obtaining this structure, and none of this 
literature (to the authors' knowledge) addresses a general methodology for obtaining a minimal M-A model for a wide 
class of uncertainty. Since having a A matrix of minimum order would improve the efficiency of structured singular 
value (or multivariable stability margin) computations, a method of obtaining a minimal M-A model would be useful. 
This paper presents a generalized method of obtaining a minimal M-A structure for systems with real parameter 
variations. 


1. Introduction 

Robust control theory for both analysis and design has been the subject of a vast amount of research in this 
decade [1-35]. In particular, robust stability and performance have been emphasized in much of this work, as, for 
example, in the development of H°° control theory [10-15, 19-23]. Moreover, the development of robust control 
system design and analysis techniques for unstructured [1-9, 13, 19, 21] as well as structured [16-35] plant uncertainty 
continues to be the subject of much research - particularly the latter. Unstructured plant uncertainty arises from 
unmodeled dynamics and other neglected phenomena, and is complex in form. This uncertainty is called "unstructured" 
because it is represented as a norm-bounded perturbation with no particular assumed structure. Plant uncertainty is 
called "structured" when there is real parameter uncertainty in the plant model, or when there is unstructured complex 
uncertainty occurring in the system at multiple points simultaneously. Plant parameter uncertainty can arise from 
modeling errors (which usually result from assumptions and simplifications made during the modeling process and/or 
from the unavailability of dynamic data on which the model is based), or from parameter variations that occur during 
system operation. 

Robust control design and analysis methods for systems with unstructured uncertainty is accomplished via 
singular value techniques [1-9, 21 ]. For systems with structured plant uncertainty, however, the structured singular 
value (SSV) [16-27 ] or multivariable stability margin (MSM) [28-33 ] must be used. In order to compute the 
SS V or MSM, the system is usually represented in terms of an M-A model. The "M" represents the transfer function 
matrix M(s) of the nominal system, and "A" represents an uncertainty matrix acting on M(s). In general, A is a block 
diagonal matrix, and for real parameter uncertainites the diagonal elements are real. As indicated in the literature 
[17,18,20,28 ], this structure can always be formed for any linear interconnection of inputs, outputs, transfer 
functions, parameter variations, and perturbations. However, very little of the literature discusses methods for 
obtaining an M-A model. For unstructured uncertainites, this model is very easy to obtain. However, for real 
parameter variations, forming an M-A model can be very difficult. In [29 ], De Gaston and Safonov present an M-A 
model for a third-order transfer function with uncertainty in the location of its two real poles and in its gain factor. 
Although the given M-A model is easily obtained for this simple example, other examples do not yield such a 
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straight-forward result. A general state model of M(s) for additive real perturbations in the system A matrix (where A 
is assumed to be closed-loop) is discussed in [34], Unfortunately, this model is not general enough for many 
examples, since system uncertainty is restricted to the A matrix and the uncertainty class is restricted to be linear. 
Morton and McAfoos [26 ] present a general method for obtaining an M-A model for linear (affine) real perturbations 
in the system matrices (A,B,C,D) of the open-loop plant state model. In this model, an interconnection matrix P(s) is 
constucted first for separating the uncertainties from the nominal plant, and then M(s) is formed by closing the 
feedback loop. The M-A model thus formed can be used in performing robustness analysis of a previously determined 
control system. If the feedback loop is not closed, p-synthesis techniques [19-21] can be applied to the M-A model for 
robust control system design. Morton's result essentially reduces to that of [34] when the perturbations occur only in 
the A matrix (and the A matrix of [34] is assumed to be open-loop). An algorithm for easily computing M(s) based 
on Morton's result is presented in [35]. Although this method of constructing an M-A model is general for linear 
uncertainties, many realistic problems require a more general class of uncertainties. Furthermore, no consideration is 
given to obtaining a minimal M-A model, where "minimal" refers to the dimension of the A (or M) matrix. Since the 
M-A model is a nonunique representation, it would be prudent to obtain one of minimal dimension so that the 
complexity of the SSV or MSM computations during robust control system design or analysis could be minimized. 
However, none of the literature (to the authors' knowledge) addresses the issue of minimality. 

This paper presents a methodology for constructing a minimal M-A model for systems with real parametric 
multilinear uncertainties, where the term "multilinear" is defined as follows: 

Definition : A function is multilinear if the functional form is linear (affine) when any variable is allowed to vary 
while the others remain fixed. For example, f(a,b,c) = a + ab + be + abc is a multilinear 
function. 

Thus, the allowance of multilinear functions of the uncertain parameters provides a means of handling cross-terms in 
the transfer function coefficients. A procedure is proposed for obtaining this model in state-space form for uncertain 
single-input single-output (SISO) systems, given the system transfer function in terms of the uncertain parameters. 

An extension of this result to multiple-input multiple-output (MIMO) systems will be given in a subsequent 
publication. In this development, M(s) will represent the nominal open-loop plant, so that the resulting M-A model 
may be used for robust control system analysis or design. The state-space form used in modeling M(s) is an extension 
of Morton’s result for real parametric linear (affine) uncertainties [26]. The paper is divided into the following 
sections. A formal statement of the problem to be solved in this paper is presented in Section 2, followed by a 
discussion of minimality considerations in Section 3. The approach is presented in Section 4, a proposed solution to 
the problem is presented in Section 5, and the proposed procedure for finding a minimal M-A model is summarized in 
Section 6. Several examples demonstraing the proposed solution are given next in Section 7, followed by some 
concluding remarks in Section 8. 


2. Problem Statement 

(liven the transfer function of an uncertain system, G(s,8), in either factored or unfactored form, as a 
function of the uncertain real parameters, 6 , find a minimal M-A model of the form depicted below in Figure 1: 


G(s, 8) 



Figure 1. Block diagram of the General M-A Model 


such that: 


1. The diagonal uncertainty matrix. A, is of minimal dimension. 

2. The model of the nominal plant, M(s), is in state-space form. 


356 



The model must handle multilinear uncertainty functions in any or all of the transfer function coefficients. In 
order to construct a minimal M-A model, the dimension of the A matrix must be minimized. Hence, factors which 
have been found to affect the dimension of the M-A model will be discussed next, followed by the approach used in 
forming a solution to this problem. 


3. Minimality Considerations 

In constructing an M-A model of an uncertain system, the A matrix can become unnecessarily large due to 
repeated uncertain parameters on its main diagonal. It is therefore of interest to examine the factors which can cause 
this repetition, so that the number of repeated uncertain parameters can be minimized. A factor which can be shown to 
cause unnecessary repetition in the A matrix is the particular realization used in representing the system. Examples 
can be constructed which demonstrate this effect, and it appears that a cascade realization (and, in particular, cascaded 
uncertain real poles and zeros) is a desirable form for obtaining a minimal M-A model. Thus, a general cascade-form 
realization will be part of the approach taken in constructing a minimal M-A model. A problem arises, however, in 
that some transfer functions have a form which precludes cascading uncertain real poles or zeros, such as: 

G(s,S) = b l s2 ± b? s + b 3 G(s,S) _ (s + 8i)(s + 82) 

(s + 0i)(s + 02) f s 2 + ai s + a2 

where 0 * and ©2 are assumed here to be uncertain (and hence a function of 8). Cascading the poles and zeros for either 
case would result in improper transfer function blocks to be realized. For these cases, it is unavoidable for the 
minimal A matrix to have repeated uncertain parameters on the main diagonal. However, for each inseparable pole or 
zero pair it is only necessary to repeat one uncertain parameter. This issue will be addressed in the proposed solution, 
and a minimal M-A model for the First transfer function above will be given as an example. 

Another factor which affects the dimension of the M-A model is the form of the coefficients in the system 
transfer function. If any of the coefficients is a nonlinear function of the uncertain parameters instead of a multilinear 
function (e.g., there are squared uncertain terms in any of the coefficients), then extra dependent uncertain parameters 
must be defined in order to represent these terms in a multilinear form. For example, would be represented as 
bi&2 where 82 = 8j, and both 8j and 82 would appear in the A matrix. Thus, for this case, it is again necessary that 
the minimal A matrix contain repeated uncertain parameters on its main diagonal. An example illustrating this 
situation will be presented later. 

These issues are addressed in the proposed solution for constructing a minimal M-A model. The approach 
taken in forming this solution is described in the next section. 


4. Approach 

Based on the problem definition and the minimality considerations outlined above, several issues will be 
addressed in forming a solution to the problem of constructing a minimal M-A model given the transfer function of an 
uncertain system. First, a general cascade-form realization will be found which can be used to obtain a minimal M-A 
model. Second, the minimal A matrix will be determined for any uncertain system such that extra dependent 
parameters are assigned to account for inseparable pairs of uncertain real poles or zeros as well as non-multilinear (e.g., 
squared) terms. Third, a method of obtaining a state-space realization of M(s) for any uncertain system will be found. 
Therefore, the proposed approach for constructing a minimal M-A model is given as follows: 

1 . Obtain a cascade-form realization of the system so that the state-space uncertain model can be written as: 

x = A x + B u (1) 

y = C x + D u 

where: A = A q + [AA], B = B q + [AB], C = C q + [AC], D = D q + [AD] ( 2 ) 

The terms with the "o" subscript (A q , B q , C , D ) represent the nominal matrix components, and the "A" terms 

(AA, AB, AC, AD) represent the uncertain matrix components. To eliminate confusion of the A notation, the 
diagonal uncertainty matrix. A, of the M-A model will be represented as [A], and the AA, AB, AC, and AD matrices 
will be represented as [AA], [AB], [AC], and [AD] wherever clarification is required. 
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2. Obtain a minimal M-A model as described in the problem definition and pictured in Figure 1, where: 
a. The minimal uncertainty matrix, A, is defined as: 

A = diag [ 8j, 5 2> § m ] = diaglfipfip] = diag[5] (3) 

where: A € R mxm , 8i e R m| , 8 D e R mD , 6 e R m 

and: m = the minimal number of uncertain parameters 

mj = the number of independent parameters given in G(s,8) 
mj) = the minimal number of dependent (or repeated) parameters 

Also: p = [A] q (4) 

P r* DlTlp 

q = the uncertain variables output from M(s), Q € 

Since an M-A model is minimal if the dimension of the A matrix, m, is minimal, where m depends 
on mj and m D , with mj being given and fixed, a formal definition of a minimal M-A model can be 

stated as follows: 


Definition : An M-A model is minimal if m D - i.e., the number of dependent (or repeated) parameters in the A matrix - 
is minimal (or zero, if possible). 


b. The state-space model of the nominal plant, M(s), is an extension of Morton's result [26] and has 
the following form: 


x 


Ao x + [ B,p B„ ] [ J 



( 5 ) 


where B X p, Cq X , Dqp, Dq U , and Dyp are constant matrices. Thus, M(s) can also be written in the 
equivalent shorthand notation defined as follows: 



M ll(s) 

m 12 (s) 


A o 

B xp 

B o 

M(s) = 

M 21 (s) 

m 22 ( s ) 


C 

L c o 

D <p 

D yp 

D qu 

°o _ 


where: Mj i(s) 

M 12 (s) 
M 2 i(s) 
M 22 (s) 


Q( s ) / P( s ) = Cq X (si - Aq) 1 B X p + Dqp 
y(s) / p(s) = Cq X (si - Aq) ^ Bq + Dq U 
<l(s) / u(s) = C 0 (si - Aq) ^ B X p + Dyp 
y(s) / u(s) = C 0 (si - Aq )' 1 Bq + D 0 


( 6 ) 


( 7 ) 


It should be noted that in [26] the D^ matrix was required to be zero. In this paper, however, D^ is 
allowed to be nonzero in order to model the multilinear (cross-product) uncertain terms. 

The results for constructing a minimal M-A model via this approach are presented in the next section. 
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5. Proposed Solution 


The proposed solution will be presented in two parts: the results for obtaining a cascade-form realization of 
the uncertain system will be summarized first, followed by the results for obtaining the state-space realization of a 
minimal M-A model. 

5.1 Cascade-Form Realization 

Given the transfer function of an uncertain system in terms of its uncertain parameters, G(s,8), it is desired to 
realize the system in a cascade form of First- and second-order subsystems. Thus, if the transfer function is given in 
unfactored form, the numerator and denominator polynomials must be factored into first- and second-order terms. The 
given transfer function will then be represented as follows: 

G(s,8) = K y (8) G c (s,S) G r (s,8) K u (8) (8) 

where K u and K y represent input and output gain terms, respectively, and Gr and G^ represent the real (first-order) and 
complex (second-order) transfer function components, respectively. Then : 


Gr(s,8) = GR k (s,8) GR k j(s,8) ... Gr 2 (s,5) Gr^s.S) 

(9) 

Gc(s,8) = Gcj(s,8) Gq ^(s.S) ... Gc 2 (s,8)Gcj(s,S) 

(10) 

Gr. ( S ,8) = P 2M s + P 2i 


S + cti 

(11) 

Gc (S 8) — t>3i — 2 S 2 + b3i-l S + b3i 


s 2 + a 2 i i s + a 2 i 

(12) 


and: k = number of real (first-order) blocks 

1 - number of complex (second-order) blocks. 

Any or all of these transfer function coefficients may be uncertain. The uncertainty may arise from either the 
coefficient itself being uncertain, or from the coefficient being a function of one or more uncertain variables. 
Therefore, for either case, any of the coefficients may be a function of 8. Furthermore, the uncertain variables may 
have either an additive or multiplicative form: 

e - e Q + ^ , e = e Q (1 + Sfe) (13) 

The following cascade-form state-space realization of this system is proposed: 



1 

> 

50 

o 

b rk u 

o 

W' 

II 

b c c r ac 

EC ER K u 


Ky D c C R K y c c 

K y c C d R K U 


where: 
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a r = 


A Rl 

Br 2 Cri 
b r D r c R 

R 3 k 2 1 

B R D d D h Cd 
K 4 r 3 R 2 R 1 


B R, C r 2 


Rc. r 2 


BR k .! DR k-2- - D R 2 C Rl V 1 Dr m" ° R 3 Cr 2 ' 
BR k DR k-r DR 2 CR i B Rk DR k-r DR 3 CR 2 ' 


B R k CR k-l A R k. 


B R = 


b R 3 Dr 2 D Ri 

Br „ Dr 3 Dr 2 D R, 


BR k-l Dr Ic-2 ' Er 2 D r , 

_ BR k DR k.l- DR 2 DR l _ 

Cr ■ [ DR k DR k.l" DR 2 CR l DR k DR k-r- D R 3 C R 2 ■■■ DR k CR k-l C Kk] (17) 

Dr * |^ DR k DR k .i ' Dr 2 d Ri1 (18) 

The A q . B(\ C(\ and matrices have the exact same form as (15) - (18), except that the subscripts "R" and ”k” are 

replaced by "C” and ”1", respectively. The submatrices are defined as follows: 

A Ri = -<Xi B Ri = 1 

(19) 

C Ri = p2i - Oi foi l °Ri = foi l 

-L\ Bo = [?l 

(20) 

Cq = [ (foi - »2i foi-2) (foi l - a 2i-l foi-2) ] D G = foi-2 

The realizations { A R ., B R . , C R ., Dr. ) and { Aq, Bq, Cq, Dq ) represent the 1 th real (first-order) and complex 
(second-order) systems Gr.(s,8) and Gq(s,8), respectively. Thus, for the real subsystems, i = 1, 2, , k, and for the 
complex subsystems, i = 1, 2 , .... 1. 
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The resulting cascade-form realization of the uncertain system is therefore given from (14) as: 

BrK u 
B e Dr K u . 

( 21 ) 

[ y Dc Dr K u 

The above model is a general cascade-form realization for any uncertain open-loop SISO transfer function. 

The model does not, however, handle nonmonic denominator polynomials with uncertain leading coefficients. This 
would result in fractional (i.e., rational) matrix elements in the realization with uncertain parameters in the 
denominator of these elements. For real uncertain poles or zeros, two factors determine whether the real (first-order) or 
complex (second-order) block form should be used The first is the nature of the uncertainty associated with these 
terms, and the second is the form of the transfer function. If the real pole or zero locations are the uncertain parameters 
and the transfer function form allows these poles or zeros to be separated out, then the real block form should be used. 
If the transfer function form does not allow this separation, then the complex block form must be used. Furthermore, 
if there is a pair of uncertain poles or zeros that cannot be cascaded, then the resulting minimum A matrix will have a 
repeating parameter on the main diagonal for each inseparable pole or zero pair. Alternatively, if the coefficients of the 
second-order polynomial associated with the real poles are the uncertain parameters, then the complex block form 
should be used. These cases will be illustrated in the Examples section of this paper. The formulation of the 
minimal M-A model will be presented next. 

5.2 Mini m al M -A Mo del 

In formulating the minimal M-A model, the minimal A matrix must be determined first, followed by the 
state-space realization of M(s). Thus, the results for formulating this model will be presented in this order. 

5.2.1 Minimal A Matrix 

The minimal A matrix is defined as in (3) with: 

m = mj + m D (22) 

where mj is the number of independent uncertain parameters, and mj) is the number of dependent uncertain parameters 
that must be added. The uncertain independent parameters are those defined in G(s,8). However, as discussed 
previously, the dependent uncertain parameters are those independent parameters that must be repeated due to non- 
multilinear terms in the transfer function coefficients and/or pairs of uncertain real poles or zeros that cannot be 
cascaded. Thus, for A to be minimal, m D (or 8 D ) should be minimized. It can be shown that if the system transfer 
function is formed from a given minimal M-A model of an uncertain system, the coefficients of the numerator and 
denominator polynomials will be multilinear functions of the uncertain parameters. Unfortunately, the converse is not 
necessarily true in general because of the dependence of the M-A model on the realization used for the plant. If the 
general cascade-form realization posed in this paper is used, however, the multilinear form of the transfer function 
coefficients can be used to establish that m = mj (i.e., m D = 0), unless there are real uncertain pairs of poles or zeros 
that cannot be cascaded. Furthermore, it can be shown that if the coefficients of all the factors of the numerator and 
denominator polynomials are multilinear functions, then the coefficients of the expanded polynomials will also be 
multilinear. However, if there are non-multilinear uncertain terms in the transfer function, then dependent parameters 
must be defined (and added to A) to represent the non-multilinear term in a multilinear form. Moreover, if the non- 
multilinear term is of the form 8 11 , then n-1 dependent parameters must be defined. If there are pairs of real uncertain 
poles or zeros that cannot be cascaded, then one additional dependent parameter must be added for each pair, and the 
dependent parameter can be either of the uncertain real parameters in the pair. Therefore, the number m, as determined 
by these rules, is the minimal dimension of the A matrix for the uncertainty class considered in this paper. Once this 
minimal dimension is determined, the A matrix can be defined as a diagonal matrix, as in (3), with the specified 
uncertain parameters on the main diagonal. Examples which illustrate these cases will be presented later in Section 6. 

5.2.2 State-Space Realization of M(s) 

Once the cascade -form realization has been determined, the system can be modeled as in (1) and (2), where 
[AA], [AB], [AC], and [AD] are known functions of the uncertain parameters. Since any non-multilinear terms have 


A = I 0 A * 

L Be Cr 

c = [ Ky D c Cr 


0 

Ac 
KyCc ] 


B = 
D = 


361 



been redefined in a multilinear form when the minimal A matrix is determined, these matrices are multilinear functions 
of the parameters. In order to obtain a state-space model for M(s) as defined in (5), expressions for these uncertainty 
matrices must be determined in terms of the matrices B X p, Cq X , Dqp, Dq U , and Dyp from the model. Using (4) and 
(5), these expressions are determined as follows: 

[AA] = B xp [A ]( I - Dqp [A J )’ 1 C qx = B xp ( I - [A ] Dqp Y l [A ] C qx 

[AB] = B X p [A ] ( I - Dqp [A ] ) ^ D qu = B X p ( I - [A ] Dqp ) ^ [A ] D qu (23) 

[AC] = Dyp [A ] ( I - Dqp [A ] ) _1 C qx = D yp ( I - [A ) Dqp Y l [A ] C qx 

[AD] = D^ [A] (I- Dqp [A ])' 1 D qu = Dyp ( I - [A ] Dqp Y l (A ] D qu 

The inverse term makes computation of Dqp very difficult. Furthermore, the matrix inversion can cause [AA], [AB], 

[AC], and [AD] to have fractional (i.e., rational) elements with uncertain parameters in the denominator, which is not 
allowed in the uncertainty class being considered. Thus, it is desirable to represent this term in expanded form as 
follows: 


( I - [A ] Dqp )-! = I + [A ] Dqp + ( [A ] Dqp ) 2 + ( [A ] Dqp ) 3 + . . . (24) 

where the latter form in (23) has been assumed. Then the above equations can be rewritten as: 

[AA] = B X p[A]Cq X + B X p{[A]Dqp + ([A]Dqp ) 2 + ([A]Dqp ) 3 + ... } [A ] C qx 

[AB] = B X p [A ] Dq U + B X p { [A ] Dqp + ( [A ] Dqp ) 2 + ( [A ] Dqp ) 3 + ... } [A ] D qu (25) 

[AC] = Dyp [A ] C qx + Dyp { [A ] Dqp + ( [A ] Dqp ) 2 + ( [A ] Dqp ) 3 + . . . } [A ] C qx 

[AD] = Dyp [A ] Dq U + Dyp { [A ] Dqp + ( [A ] Dqp ) 2 + ( [A ] Dqp ) 3 + ... } [A ] D qu 

The second group of terms add in the cross-terms of the multilinear uncertainty functions. Each term in the series adds 
a higher-order cross-product term. Since [AA], [AB], [AC], and [AD] are multilinear functions with a finite number of 
terms, the Dqp matrix can be defined to have a special nilpotent structure such that: 

Dqp r+1 =0 (26) 

and: ( I - [A ] D qp Y l = I + [A ] Dqp + ( [A ] D qp ) 2 + . . . ( [A ] Dqp ) r (27) 

where r is the order of the highest cross-term occuring in [AA], [AB], [AC], and [AD], i.e.: 

r = max ( Oa, Or, Oc, Oj) ) (28) 

where O^, Og, O^, and Op represent the order of the highest-order cross-product term in [AA], [AB], [AC], and [AD], 
respectively. Cross-product term order is defined as: 

order ( 8 j 82 S 3 . . . 8 j ) = i - 1 (29) 

where i = 1, 2, . . . , m. Thus, the maximum value of r is r max = m-1, where m is the dimension of the A matrix. 
The required structure for Dqp to satisfy (26) and (27) is given as follows: 

1 . ) dii = 0 ; i = 1 , 2 , . . . , m 

2. ) If djj * 0 , then for i = 1 , 2 , . . . , m and j = 1, 2 m: (30) 

a. ) dj, = 0 ; 

b. ) d i©lj©l ■ 0 or d i® 2 j ®2 = 0 or ... or dje^.^j^^.i) = 0 
where the symbol ”©" represents "modulo m” addition. The desired equations can therefore be written as: 
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(31) 


[AA] — B X p [A ] Cq X + B X p { [A ] Dqp + ( [A ] Dqp )^ + . . . ( [A ] Dqp ) r } [A ] Cq X 

[AB] = B xp [A ] D qu + B xp { [A ] Dqp + ( [A ] ) 2 + . . . ( [A ] Dqp ) r ) [A ] D qu 

[AC] = D yp [A ] C qx + D yp { [A ] Dqp + ( [A ] Dqp )2 + . . . ( [A ] Dqp >r } [A ] C qx 

[AD] = Dyp [A ] D qu + Dyp { [A ] Dqp + ( [A ] Dqp ) 2 + ... ( [A ] Dqp f } [A ] D qu 

where "r" is defined in (28). Since the [AA], [AB], [AC], and [AD] matrices are known for the given system, the 
equations in (31) are used to determine B X p, Cq X , Dq U , Dyp, and Dqp. Once these matrices are obtained, the state- 
space model of M(s) is found. Hence, a minimal M-A model has been formed. 

A procedure which summarizes the necessary steps in obtaining a minimal M-A model using these results is 
presented next 


6. Summary of.. Procedure 


The following is a summary of the procedure implied by the above proposed approach for forming a minimal 
M-A model of a given uncertain system: 

i. ) Obtain the system transfer function in factored form. The coefficients of each factor should be a multilinear 

function of the uncertain parameters. If necessary, define new dependent parameters to represent any non- 
multilinear terms in a multilinear form. 

ii. ) Define the number of parameters in the A matrix, m, using (22). In so doing, determine if any new parameters 

are required to model inseparable uncertain real pole or zero pairs. If there are inseparable real pairs, either 
uncertain parameter in the pair may be repeated. 

iii. ) Define the minimal A matrix as in (3), using the independent parameters defined in the given transfer function 

as well as those defined in steps i.) and ii.) above. 

iv. ) Obtain a cascade-form realization for the system as a function of the uncertain parameters. 

v. ) Express the system matrices as in (2). 

vi. ) Determine the maximum order of cross-product terms, r, in [AA], [AB], [AC], and [AD] as defined by (28) and 

(29). Then [AA], [AB], [AC], and [AD] have the form represented in (31), where Dqp has the special 
(nilpotent) structure summarized by (30). 

vii. ) Express [AA], [AB], [AC], and [AD] as: 

[AA] = [AA 0 ] + [AA^ + [AA 2 ] + . . . + [AA r ] 

[AB] = [AB 0 ] + [AB!] + [AB 2 ] + ... + [AB r ] (32) 

[AC] = [AC 0 ] + [AC i] + [AC 2 ] + ... + [AC r ] 

[AD] = [AD 0 ] + [ADj] + [AD 2 ] + . . . + [AD r ] 

where the subscript i represents the cross-terms of order in each uncertainty matrix. 

viii. ) The B X p, Cq X , Dyp, and Dq U matrices are found using the expansion described in [26] for the uncertainty 

matrices having zero-order cross-product terms; i.e. define: 


AAo AB 0 
_AC 0 AD 0 


= Mi 5i + M 2 f>2 + - • • + M m 5m 


(33) 


where the Mj matrices are appropriately partitioned. For the case of repeated parameters due to inseparable real 
poles or zeros, the Mj matrix associated with the repeated parameter must be nonzero . These matrices can be 
decomposed into the product of appropriately partitioned column and row matrices as follows: 
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Mi = 


M Bi 

Mdu 


[MqM^] 


(34) 


where Mg. forms the i* column of B xp , M D] . forms the i* column of Dyp, Me- forms the i th column of 


^ orms ^ column of Dq U . Thus: 


Bxp — [ 

Mg, 

Mg 2 

Mr.] 

Dyp = [ 

MDj! 

m D i2 • • • 

Md,J 

II 

cT 

Mq T 

Me/ ••• 

Mc„] T 

< 

11 

MD2, 7 

MdJ ••• 

Mrvl 7 


(35) 


ix. ) Use the higher-order cross-terms of [AA], [AB], [AC], and [AD], as in (32), to determine the elements of the 

D qp matrix. Begin with the first-order terms and specify as many elements as possible. Continue with the 
second-order terms, and proceed until all elements of Dqp are specified. Check Dqp to ensure that the required 
special structure of (30) and, hence, (26) is satisfied. 

x. ) Form the minimal M-A model as given in (3), (5) and (6), and as pictured in Figure 1. 

It should be noted that the matrices Mg., Mq, M^., and obtained in decomposing the Mj matrices 

in (34), are not necessarily unique. A method of formalizing this decomposition for computer implementation will 
not be addressed in this paper. However, an algorithm is presented in [35] which accomplishes this decomposition for 
a more restrictive uncertainty class. Some examples will be given next to illustrate these results. 


7. Examples 

The following examples illustrate the proposed procedure presented above. Due to space limitations, details 
of each step will not be included. However, the results that are presented for each example should be fairly easily 
obtained. 

Example 7.1 

Consider the following uncertain system: 

(Pi S + P 2 XP 3 S + p 4 )(bl S 2 + t >2 S + t> 3 ) 


G(s,8) = 


(s + aO(s + ot 2 )(s 2 + ai s + a 2 ) 


where: 04 = a lo + 5 ai , cc 2 = ot2 0 + 8 a 2 • a l = a l 0 + 8a l ’ 32 = a 2 0 + 8a 2 

Pl=Pl 0 + 8 Pl • P2 = P2 0 + 8 P 2 • P3 = P3 0 + 5 P 3 • P4 = P4 0 + 8 p 4 

b l = b l 0 + 6b l ’ b2 = b 2 0 + 5 b 2 • b 3 = b 3 0 + 8 b 3 

The cascade-form realization of this example is found in a straight-forward manner to be: 


A = 


-ai 0 0 0 

B = 

1 

Pi 

p 2 -aiPi -<X2 0 0 

0 

0 0 0 1: 


. p3(P2-«lPl) P4-a2P3 * a 2 * a l 


1 

5 s 

XP 

u> 

1 
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Example 7.2 

This example illustrates the case of an inseparable uncertain real pole pair. Consider the following system: 

G(s, 8) = b i s2 + ^ s ±- k 

(s + 0i)(s + 02) 

where: bj = bi Q + 85^ , b^ = b2 0 + fy>2 * b 3 = b 3 G + ^b3 

6l = e l o + 5 0l • °2 = 0 2 o + 5 02 


Since the numerator is second order with uncertain coefficients, the uncertain real poles in the denominator cannot be 
separated into the real cascade form. The denominator must therefore be expanded, and the complex (second-order) 
block used in the realization, which is given as follows: 


A = 


0 1 
-0102 *(01 + 02 ) 


, B = 

0 


. 1 . 


c = [ (b 3 - 0!0 2 bi) (b 2 - (0i+0 2 )bi) ] . D = M 


Since there is one inseparable uncertain pair of poles, either Sq ^ or Sgj must be repeated in the A matrix. (It can be 
shown that if this is not done, Dqp will not have the required structure and hence the higher-order cross-product terms 
will not be modeled correctly.) Since there are no non-multilinear terms in any of the transfer function coefficients, 
m = 6 (i.e., five given independent uncertain parameters plus one dependent parameter for the inseparable pole pair). 
The resulting A matrix can therefore be defined as follows: 

A = diag [8e lt 8b 3 l 


where 8 q 1 was 


arbitrarily chosen to be repeated. Using the proposed procedure, the following results are obtained: 


Bxp 


0 0 0 0 0 O' 

-1 -1 -1 0 0 0 . 


Dyp — [-bio -bio -bio 1 1] 



0 2 O 

0 

010 

01o02o 


0 

1 


0 

1 

1 

(01o+02o) 

1 

0 


, Dqu — 


0 

0 

0 

-1 

0 

0 - 


Dq) - 


0 

0 

l/02o 

1 

0 

0 


0 

0 

0 

1 

0 

0 


0 0 

0 0 

0 0 

1 0 

0 0 

0 0 


0 0 
0 0 
0 0 
0 0 
0 0 
0 0 
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E x ample 7, 3 


This example illustrates the case of non-multilinear terms in the transfer function. Consider the second-order 
uncertain system: 


G(s,5) = — i - 

s 2 + 2a s + (o 2 , where: a = a+ 8 , (o = co + 5 

o a o co 

This is a second-order system with uncertain complex poles. The uncertainty appears in the real and imaginary 
components of the complex poles. The constant coefficient, co 2 , is not a multilinear function of the uncertain 
parameters. Substituting for a and co in the above transfer function yields the following equation: 

s 2 y = u - 2a sy - co 2 y - 25 sy - ( 2co 5 + 5 2 ) y 
o o 3 a 7 ' o co co 

In forming an M-A model, the problematic term is 5^ 2 because it is not multilinear. In order to represent this 
equation in multilinear form, the following dependent variable is defined: 


\ = 5 

3 co 


so that: 


- u ■ 2a, sy ■ o> 0 2 y • 28, sy - (Zo^+S^y 


Then the following equations can be defined: 
q l = -2 sy 

q 2 = -y 


P 1 - S c q l 
p 2 - 


p 3 " S 3 q 3 


q 3 = - 20 ) o y + p 2 
Thus: s^y = u - 2a q sy - co^ y + 

The realization of M(s) for the resulting M-A model can be depicted as follows: 


Xl 



0 1 


Xf 

+ 

f° 0 01 

Pi' 

P2 

+ 

0 

X2. 


. -CC& -2 CTo . 


-X2. 


1 0 1 J 

_P3_ 


. 1 . 


-qr 


' 0 -2 ' 


Y, 1 

0 0 O' 


Pl' 

Q2 

= 

-1 0 


Xl + 
Y ^ 

0 0 0 


P2 

_<13_ 


. -2co 0 0 


L *2 J 

-0 1 0- 


.P3_ 



The A matrix is given by A = diag [8 ,8 , 8, ], where 8„ = 8 . 

a to 3 3 to 


These examples illustrate the proposed procedure for forming a minimal M-A model of an uncertain system. 
Although all the steps involved in obtaining these results have not been included, the stated results should provide a 
guide in performing the steps of the proposed procedure. It should be noted that, for ease of hand computation, the 
examples included only the simplistic (and less realistic) case in which the coefficients themselves are the uncertain 
parameters. However, it is emphasized that the proposed procedure does handle the more realistic case in which the 
uncertain transfer function coefficients are multilinear functions of the uncertain parameters. 
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8. Concluding Remarks 


This paper has presented a proposed procedure for forming a minimal M-A model of an uncertain system 
given its transfer function in terms of the uncertain parameters. The uncertainty class considered in this paper allows 
the transfer function coefficients to be multilinear functions of the uncertain parameters, and the uncertainties may arise 
in the A, B, C, and D matrices of the system model. The proposed procedure involves realizing the system in a 
cascade form, determining the minimal A matrix of uncertain parameters, and obtaining a state-space model for the 
nominal system, M(s). Three examples were given to illustrate the proposed procedure. The first example had eleven 
independent uncertain parameters, which arose in the A, B, C, and D matrices of the system realization. The second 
example had uncertain parameters arising in the A, C, and D matrices only. This example illustrated the formulation 
of a minimal M-A model for a system with inseparable real uncertain poles, and involved repeating an uncertain 
parameter in the A matrix. The last example had uncertainty in the A matrix only, and illustrated a method for 
handling non-multilinear terms. 

Further work on the proposed procedure will be to include systems having a nonmonic characteristic 
polynomial with an uncertain leading coefficient, as well as systems having an inner feedback loop which may or may 
not have uncertainties. The latter case may require a modification in the formulation of the cascade realization. 
Although the procedure presented in this paper is for SISO systems, an extension to MIMO systems will be 
forthcoming, and should primarily involve modifying the cascade-form realization. Other areas of future work include 
development of a simple means of verifying the minimality of a given M-A model, and development of a method of 
reducing a nonminimal M-A model to a minimal form. 

REFERENCE LIST 


tl] 

[2.1 

[3.1 

[4.] 

[5.] 

[ 6 .] 

[7.1 

[8.1 

[9.1 

[ 10 .] 

[ 11.1 

[ 12 .] 

[13.] 

[14.1 

[15.1 

[16.1 


J. C. Doyle and G. Stein, "Multivariable Feedback Design: Concepts for a Classical/Modem 
Synthesis", IEEE Trans, on AuL Contr., Vol. AC-26, pp.4-16, Feb. 1981. 

M. J. Chen and C. A. Desoer, "Necessary and Sufficient Condition for Robust Stability of Linear 
Distributed Feedback Systems", Int. J. Control, Vol. 35, No. 2, pp. 255-267, 1982. 

N. A. Lehtomaki, N. R. Sandell, and M. Athans, "Robustness Results in Linear-Quadratic-Gaussian Based 
Multivariable Control Design", IEEE Trans, on Aut. Contr., Vol. AC-26, pp. 75-92, Feb. 1981. 

N. A. Lehtomaki, D. A. Castanon, B.C. Levy, G. Stein, N. R. Sandell Jr., and M. Athans, 

"Robustness and Modeling Error Characterization", IEEE Trans, on Aut. Contr., Vol. AC-29, No. 3, 

March 1984. T _ 

I. Postelthwaite, J. M. Edmunds, and A. G. J. MacFarlane, "Principal Gains and Principal Phases In The 
Analysis of Linear Multivariable Feedback Systems", IEEE Trans, on Aut. Contr., Vol. AC-26, pp. 32-46, 
1981. 

M. G. Safonov, A. J. Laub, and G. L. Hartmann, "Feedback Properties of Multivariable Systems: The Role 
and Use of the Return Difference Matrix", IEEE Trans, on Aut Contr., Vol. AC-26, pp. 47-65, 1981. 

J. B. Cruz, J. S. Freudenberg, and D. P. Looze, "A Relationship Between Sensitivity and Stability of 
Multivariable Feedback Systems", IEEE Trans, on Aut Contr., Vol. AC-26, pp. 66-74, 1981. 

J. C. Doyle, "Robustness of Multiloop Linear Feedback Systems", 17th IEEE Conf. on Decision and 

Control, San Diego, USA, Jan. 1979. . . . 

J. S. Freudenberg, D. P. Looze, and J. B. Cruz, "Robustness Analysis Using Singular Value Sensitivities , 

Int. J. Control, Vol. 35, pp. 95-116, 1982. 

B. A. Francis, J. W. Helton, and G. Zames, "H°°- Optimal Feedback Controllers for Linear Multivariable 
Systems", IEEE Trans, on Aut. Contr., Vol. AC-29, pp. 888-900, Oct 1984. 

B.-C. Chang and J. B. Pearson, "Optimal Disturbance Reduction in Linear Multivariable Systems", IEEE 
Trans, on Aut. Contr., Vol. AC-29, pp. 880-887, Oct 1984. 

B.-C. Chang, "A Stable State-Space Realization in the Fbrmulation of H°° Norm Computation", IEEE 
Trans. Aut. Contr., Vol. AC-32, pp. 811-815, Sept 1987. 

B.-C. Chang, S. S. Banda, and T. E. McQuade, "Fast Iterative Computation of Optimal Two-Block H°°- 
Norm", IEEE Trans, on Aut. Contr., Vol. AC-34, pp. 738-743, July 1989. 

B.-C. Chang, "Size Reduction in Four-Block H°° Formulation", Proceedings of 26th IEEE Conf. on 

Decision and Control, Dec. 1987. iti 

J. C. Doyle, K. Glover, P. Khargonekar, and B. A. Francis, "State-Space Solutions to Standard H 2 and rL, 

Control Problems", 1988 American Control Conference, Atlanta, GA, June 1988. 

M. G. Safonov and J. Doyle, "Minimizing Conservativeness of Robustness Singular Values", in 
S. G. Tzefestas (ed.), Multivariable Control. D. Reidel Publishing Company, 1984. 


368 


[17.] J. Doyle, "Analysis of Feedback Systems with Structured Uncertainties", DEE Proc. Part D, Vol. 129, No. 
6, pp. 242-250, 1982. 

[18.] J. C. Doyle, J. E. Wall, and G. Stein, "Performance and Robustness Analysis with Structured Uncertainty", 
Proceedings of the 21st IEEE Conference on Decision and Control, 1982. 

[19.] J. Doyle, "Synthesis of Robust Controllers and Filters", Proc. of 22nd IEEE Conf. on Decision and 
Control, Dec. 1983. 

[20.] J. Doyle, "Structured Uncertainty in Control System Design", Proc. of the 24th DEEE Conf. on Decision 
and Control, Dec. 1985. 

[21.] J. Doyle, Lecture Notes, ONR/Honeywell Workshop on Advances in Multivariable Control, Minneapolis, 
MN, Oct. 1984. 

[22.] M. G. Safonov, "Optimal H°° Synthesis of Robust Controllers For Systems with Structured Uncertainty", 
Proc. of 25th IEEE Conf. on Decision and Control, pp. 1822-1825, Dec. 1986. 

[23.] J. C. Doyle, "A Review of p. For Case Studies in Robust Control", IFAC 10th Triennial World Congress, 
Munich, FRG, 1987. 

[24.] M. K. H. Fan and A. L. Tits, "Characterization and Efficient Computation of the Structured Singular 
Value", IEEE Trans. Aut. Contr., Vol. AC-31, pp. 734-743, Aug. 1986. 

[25.] M. K. H. Fan and A. L. Tits, "m-Form Numerical Range and the Computation of the Structured Singular 
Value", IEEE Trans. Aut. Contr., Vol. AC-33, pp. 284-289, 1988. 

[26.] B. G. Morton and R. M. McAfoos, "A Mu-Test for Real Parameter Variations", Proc. of thel985 
American Control Conference. 

[27.] R. D. Jones, "Structured Singular Value Analysis for Real Parameter Variations", 1987 AIAA Conference on 
Guidance and Control. 

[28.] M. G. Safonov, "Stability Margins of Diagonally Perturbed Multivariable Feedback Systems", Proc. of the 
20th IEEE Conf. on Decision and Control, Dec. 1981. 

[29.] R. E. De Gaston and M. G. Safonov, "Exact Calculation of the Multiloop Stability Margin", IEEE Trans, 
on Aut. Contr., Vol. AC-33, pp. 156-171, 1988. 

[30.] R. S. S. Pena and A. Sideris, "A General Program to Compute the Multivariable Stability Margin for 
Systems with Parametric Uncertainty", Proc. of the 1988 American Control Conference. 

[31.] A. Sideris and R. S. S. Pena, "Fast Computation of the Multivariable Stability Margin for Real Interrelated 
Uncertain Parameters", Proc. of the 1988 American Control Conference. 

[32.] B.-C. Chang, O. Ekdal, H. H. Yeh, and S. S. Banda, "Compulation of the Real Structured Singular Value 

via Polytopic Polynomials", Proc. of the 1989 AIAA Guidance, Navigation and Control Conference. 

[33.] H. Bouguerra, B.-C. Chang, H. H. Yeh, and S. Banda, "A Fast Algorithm for Checking the Stability of 

Convex Combinations of Stable Polynomials", Technical Report, Dept, of Mechanical Engineering, Drexel 
University, Jan. 1989. 

[34.] G. A. Hewer, C. Kenney, and R. Klabunde, "A State Space Model of Structured Singular Values", Proc. of 
the 27th IEEE Conf. on Decision and Control, Dec. 1988. 

[35.] M. K. Manning and S. S. Banda, "Algorithm to Obtain M-A Form for Robust Control", Proc. of IEEE 
Conf. on Systems Engineering, Aug. 24-26, 1989, Dayton, Ohio. 


369 



